Staggered Chiral Perturbation Theory for Heavy-Light Mesons 

C. Aubin 

Dept. of Physics, Columbia University, New York, NY 10027 

C. Bernard 

Washington University, St. Louis, MO 63130 
Abstract 

We incorporate heavy- light mesons into staggered chiral perturbation theory (SXPT), working 
to leading order in l/mg, where mq is the heavy quark mass. At first non-trivial order in the 
chiral expansion, staggered taste violations affect the chiral logarithms for heavy- light quantities 
only through the light meson propagators in loops. There are also new analytic contributions 
coming from additional terms in the Lagrangian involving heavy-light and light mesons. Using 
this heavy-light SXPT, we perform the one-loop calculation of the B (or D) meson leptonic decay 
constant in the partially quenched and full QCD cases. In our treatment, we assume the validity 
both of the "fourth root trick" to reduce four staggered tastes to one, and of the SXPT prescription 
to represent this trick by insertions of factors of 1/4 for each sea quark loop. 

PACS numbers: 12.39.Fe,12.39.Hg, ll.30.Rd, 12. 38. Go 
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I. INTRODUCTION 



The lattice can make a major contribution to the understanding of flavor physics through 
the computation of the properties of heavy-hght mesons. Among the important quantities 
to calculate are the heavy-light decay constants, form factors, and B parameters. (See 
Ref. for recent reviews.) However, the "traditional" lattice approach, using quite 

massive light quarks, is likely to produce results that have large systematic errors coming 
from the long chiral extrapolation [3]. An alternative approach that seems very promising 
is to use staggered (Kogut-Susskind) fermions for both the light valence and sea quarks 
An exact chiral symmetry for staggered quarks at flnite lattice spacing 
allows simulations to be performed at small quark mass, deep in the chiral regime. 

Staggered quarks have "taste" symmetry, a four-fold remnant of the doubling symmetry. 
Taste symmetry is broken at O(a^) in the lattice spacing, resulting in rather large discretiza- 
tion effects. Even at the smaller lattice spacing employed in current MILC simulations with 
improved staggered quarks (the "flne" set with lattice spacing a ~ 0.09 fm), the ef- 
fects are not negligible in any quantity which is sensitive to pseudoscalar meson loops. This 
has been understood in the light meson sector using staggered chiral perturbation theory 



(SXPT) 
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13l |. We now extend this program to include heavy quarks, thus merging 



heavy quark effective theory (HQET) with SXPT. We then calculate the leptonic decay 
constants, Jb and Jb^, at one loop in SXPT in both the partially quenched and full QCD 
cases. This is one of the simpler quantities to calculate both in lattice simulations 3, 3| and 
in the SXPT formalism. 

Simulations with staggered fermions such as Refs. jj, S, 0, Q, , , 13 , use the "fourth root 
trick" Q], designed to reduce the number of tastes per flavor from 4 to 1 in the continuum 
limit. The validity of the fourth root trick has not been proven, although various recent 
studies have, in our opinion, made it rather plausible [13j. In the following, we assume that 
this trick is valid, and that it can be represented in the chiral theory by insertions of factors 
of 1/4 for each sea quark loop 12, Q- employ a quark- flow analysis [3| to locate 



empi 



these sea quark loops. The "replica method" jl7[ is an alternative technique for inserting 
the appropriate factors. It is theoretically cleaner than quark flow, but usually somewhat 
more complicated to implement in practice. One would start with ns staggered flelds for 
each continuum sea-quark flavor desired, and take the limit ns 1/4 at the end. The two 
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approaches yield identical results in the current case, and in all other cases that have been 
investigated to date. 

A brief discussion of the current results, as well as the application to semileptonic form 
factors, was previously presented in Ref. A detailed description of the semileptonic 



calculation is in preparation . The SXPT forms obtained have already been used for the 
chiral (and continuum) extrapolations of numerical lattice results for heavy-light semilep- 
tonic 'sj] and leptonic 9| decays. 

In this paper the heavy quark mass tuq is taken to be large compared to Aqcdj so that 
working to leading order in the heavy quark expansion {i.e., neglecting l/mg terms) is a 
reasonable first approximation. We also do not include, in our effective theory, discretiza- 
tion errors coming from the heavy quark. We assume that such errors can be adequately 
estimated independently, using HQET as the effective theory description of the lattice heavy 



quark 



21|. However, in order to separate light and heavy quark discretization errors. 



our analysis requires that the heavy quark mass not be much greater than one in lattice 
units. The rather counterintuitive restriction arises because the heavy quark action no longer 
sufficiently suppresses doublers when amq ^ 1. Indeed, in the large amq limit, the heavy 
quark action becomes a static lattice action, where there is no suppression of spatial doublers 
at all, since a static quark's energy is independent of its momentum. When heavy quark 
doublers are insufficiently suppressed, the light and heavy quarks in a heavy-light meson 
can exchange gluons of momenta n/a, leading to "mixed" 4-quark operators (products of 
heavy and light bilinears) in the Symanzik action that violate taste. In order to avoid treat- 
ing this complication — which would inextricably mix light- and heavy-quark discretization 
errors - we do not allow ara, » 1 ■„ our analys,. 0,r the other ha,rd. typ.eal values 
of amq in practical calculations with the Fermilab [22| or NRQCD j23|] heavy quarks are 
likely to be acceptable. With this restriction, all taste- violations in the O(a^) Symanzik 
action occur solely in the light quark sector; the taste-violating terms are thus the same as 
in Refs. More details, including a discussion based on the symmetries of the lattice 

action, appear in Sec. HH 

Incorporating heavy quarks into SXPT produces a large number of terms in the chiral 
Lagrangian involving combinations of the heavy-light and the light mesons. All of these 
terms, which arise due to the taste symmetry breaking, are however non-leading: They do 
not appear in one-loop diagrams for the heavy-light quantities and enter the decay constant 
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result through a single analytic term. Non-trivial taste violations arise only in the light 
meson propagators in the one-loop diagrams. Since the taste-violating low energy constants 
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in the light meson sector have been determined from simulations and SXPT fits 
the new analytic term is the only additional fit parameter beyond those low-energy constants 
that appear in the continuum, such as the B-B*-^ coupling g^^. 

This paper is organized as follows. We first describe the Symanzik action to O(a^) in 
Sec. In] The SXPT Lagrangian, here including heavy-light mesons, is then constructed from 
a spurion analysis in Sees. Illll and IIVI where the latter focuses on the O(a^) terms involving 
heavy- light mesons. In Secl^we calculate the one-loop expression for the chiral logarithms 
that arise in the leptonic decay constant, Jb^, where x is the flavor of the light valence quark. 
We write down the final results for various cases in Sec. IVll and discuss some features of the 
low energy constants (LECs) in Sec. IVIII Finite volume effects are treated in Sec. IVlTTI We 
finish with some remarks and conclusions in Sec. IIXI 

II. SYMANZIK ACTION 

For concreteness, we consider here a heavy quark that is simulated either with NRQCD 



23[ or with the Fermilab interpretation 22] of a clover or more highly improved quark. ^ 



What is crucial for us about these possible heavy quark actions is that they describe a single 
physical fermion — doubler^ masses are assumed to be larger than the physical heavy quark 
mass by an amount of order of the cutoff, i.e., ~ 1/a. This means that the heavy quark 
fields in the Symanzik action will have no degree of freedom corresponding to "taste." 

In the limit arriQ ^ 1, however, it will not be possible to neglect the doublers with either 
the Fermilab or the NQRCD actions. That is because these actions approach a static lattice 
theory, which has no intrinsic suppression of spatial doublers: The energy of a static quark 
is independent of its three-momentum. This would not be a problem if the light quark in 



^ The formalism will also apply when the heavy quark is simulated by extrapolation up in mass from a 
lighter, conventional, lattice quark. 

^ For convenience, we use the term "doubler" to describe a heavy-quark state with one or more components 
of 3-momentum near n/a. Its "rest mass" is defined to be the energy of the state when all 3-momentum 
components are exactly ir/a or 0. Since these states are always suppressed in the cases considered here, 
we do not need to be concerned that the interpretation of such states as additional heavy particles is 
problematic: With Wilson-like or non-relativistic fermions, their energy decreases when one of the large 
components of momentum changes away from ir/a. 
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a heavy-light bound state were simulated as a Wilson-like fermion. In that case, the light- 
quark Wilson term would suppress gluon exchanges of momentum vr/a between the light and 
heavy quarks, by forcing the light quark far off shell. But emission or absorption of such a 
gluon is not suppressed by the staggered light action, and indeed is simply a taste-changing 
interaction. To keep the analysis manageable, we must assume that the heavy quark action 
suppresses these exchanges. We define Am to be the splitting between (the rest masses of) 
the physical heavy quark and the closest doublers. To be able to ignore the effect of doubler 
states in the low energy effective theory, we must have ^ ^qcd- However, for the 
Fermilab or NRQCD action. Am goes like l/(a^mQ) for aniQ ^ 1; eventually the condition 
Am ^ ^QCD will be violated as rriQ grows at fixed a. 

We wish to correct a possible misunderstanding here. When arriQ ~ 1, one might worry 
that doublers cannot be neglected relative to the physical heavy quark since both physical 
and doubler masses are formally 0{l/a). However, what is important is not (say) the ratio 
of doubler mass to physical mass, but the splitting Am between doubler masses and the 
physical mass. As long as Am ^ ^qcd, the doubler states can treated as "integrated out" 
and their effects can be summarized by higher- dimension operators in the Symanzik/HQET 
theory 0, involving only the light quarks and the physical heavy quark. Systematic 
errors coming from doublers are thus included in estimates of heavy-quark discretization 
effects from higher operators. (See further discussion below). 

It is instructive to estimate Am in some practical situations. For Wilson-type quarks, the 
rest (or "pole") mass mi of the physical quark is given at tree level by ami = lii(l + amo), 
where amo = 1/{2k) — l/{2K,crit), with k the hopping parameter and Kcru its critical value. 
For the lowest doubler, with momentum vr/a in a single lattice direction, the Wilson term 1 — 
cos(ap) becomes 2 instead of 0, and the doubler rest mass mf is given by amf = ln(3-|-amo). 
For the MILC "coarse" lattices [lO|l, a~^ ~ 1.59 GeV, hicrit ~ 0.1378, and kj,, the hopping 
parameter of a 6 quark is ~ 0.08 6 |2fil| . This gives Am = mf —nii ^ 775 MeV. Alternatively, 
using the tadpole improvement [27] parameter uq = l/{8Kcrit) gives anio = 4(Kcrjt//« — 1) and 
Am ~ 735 MeV. These values of Am, which are comparable to other masses dropped in XPT 
{e.g., nip), are probably large enough to neglect the effect of doublers in the current chiral 
analysis. However, it is clear that applying this analysis to any lattices that are significantly 
coarser than the MILC coarse lattices will be problematic. For example, at a^^ = 1 GeV, 
we estimate Am ~ 350 MeV, which is certainly not much larger than Aqcd- With Am this 
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low, any effective chiral theory becomes very complicated, and there also would be practical 

problems in separating the physical state from the doubler states in simulations. 

A similar analysis applies in NRQCD. Using a bare mass aM = 2.8 for h quarks on the 

n 

MILC coarse lattices and a stabilization parameter = 2 [6], we find Aj\/ 930 MeV, which 
should be adequate. However, with a^^ = 1 GeV, Am ^ 400 MeV, again much too low to 
omit the effects of doublers from the effective chiral theory. 

If heavy quark doubler effects may be neglected, the analysis of the Symanzik action is 



very similar to that in Ref. [28|, which considers Ginsparg- Wilson valence quarks coupled 
to staggered sea quarks. Taste violations at C(a^) appear only in four-quark operators 



composed exclusively of 
identical to those in jlll 



igh t (staggered) quarks. These terms in the Symanzik action are 
13|. "Mixed" four-quark operators consisting of the product of a 
light quark bilinear and a heavy quark bilinear do not break taste symmetry. Physically, 
operators that violate taste symmetry require momentum it /a gluon exchanges, which we 
ban from the Symanzik action when we omit doubler states from the low-energy theory. 
A more rigorous proof that follows from the symmetries of the lattice action can also be 
constructed. The proof is identical to that presented in Ref. j2s|, so we do not include it 
in detail here. The basic idea is that the continuum, physical, heavy quark fields transform 
trivially under translations by one lattice unit. Staggered bilinears with non-singlet taste, 
on the other hand, get multiplied by phase factors under single-site translations. Therefore 
mixed four quark operators cannot be translation singlets unless the light bilinear is a singlet 
under taste. 

The mixed four-quark operators are thus irrelevant from the point of view of our heavy- 
light chiral theory. Such operators are invariant under the light quark chiral symmetries 
since they do not break taste (or, trivially, flavor) symmetries and are independent of the 
light quark masses (otherwise, they would have dimension greater than 6). We therefore can 
classify the discretization errors caused by such operators as "heavy-quark errors." These 
errors can be estimated by the methods of Refs. 3, I21I. (See especially the Appendix 
of Ref. [2^ for estimates of heavy quark discretization errors in a practical case.) In fact, 
one expects that errors from mixed four-quark operators will be rather smaller than the 
more familiar ones from heavy-quark bilinear operators, since the former are 0{a'^a'^) for 
improved staggered ("Asqtad") light quarks 2^; while the latter are 0{asa) or (9(a^). It is 
important to keep in mind, however, that such error estimates are based on power counting 
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and dimensional analysis and must be considered as rough guides only. Precise quantification 
of the total discretization error will always require simulation at several lattice spacings. 

Four-quark operators built entirely out of heavy quarks are also invariant (trivially) under 
the light quark symmetries. Therefore, the arguments of the previous paragraph apply, with 
the exception that the discretization errors are now (9(a5a^): There is no suppression of the 
exchange of a single gluon with momentum vr/a between the two quarks, unlike the cases 
where one or both of the quarks has the Asqtad action. 

The relevant part of the Symanzik action may now be written as S^^ = f d^xC^^, with 

£Sk _ ^Sk ^ ^2^Sk ^ (^^^4^ ^ 

where Cf^ comes from light four-quark operators only. Cf^ is the continuum limit of the 
lattice theory with n staggered light quarks (each with 4 tastes), and a single heavy quark. 
It has the form 

n 

= 5Z 9i («f - )q,+Q{ivD)Q + £gi„ons • (2) 

Here v is the heavy quark velocity, and j runs over the light quark flavors. i2sk,4 has a heavy 
quark SU{2) spin symmetry and, in the limit that rrij ^ 0, a chiral SU^An)^ x SU{'in)R 
symmetry. 

The 0{a^) term Cf in Eq. (0) is a sum over all 4-light-quark operators in the two classes 



^6^^^^^ and S^^^^^ identified in Ref. 

Cf = Cf + Cf . (3) 

Sq consists of ten operators that do not violate continuum rotation symmetry, while the 
four operators in violate this symmetry. For n staggered fields, the operators in each 

of these classes are written down in Ref. js^. Every operator is a product of two bilinears, 
each of which has the same spin and taste. Denoting the five possible spins structures by S, 
V, T, A, P, and similarly for tastes, the operators are named by the spin ® taste of their 
bilinears: 

Cf : [S X A], [S X V], [A X S], [V X S], [P x A], [P x V^], 

[A X P], [V X P], [T X V], [T X A], [V x T], [A x T] (4) 
Cf : [T, X V,], [T, X A,l [V, x T,], [A, x T,] (5) 
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In the spin indices and taste indices are contracted separately, for example, 

[T X A] = ^Y^ qi{-f^^ O ^xb)qi qjili^f, ® ^5\)qj (6) 

where 7^^, = (1/2) [7^, 7^,], 7^5 = 7^75 (similarly for tastes, with 7^ ^^), and there are 
implicit sums over the flavor indices The operators in are more complicated: spin 
and taste indices are summed together, and two terms, with different treatment of the tensor 
indices, are subtracted to cancel taste- or spin-singlet contributions. For example. 



[T^ X A^] = ^ ^ |gi (7^^ ® ^^5) Qi qj (7^^ ® ^5^) qj - (7^^5 ® ^^5) g^g^- (75^.^ ® ^5^) qj 

(7) 

Because the same index (/i) appears four times in each term in Sq^^^\ these operators are not 
invariant under continuum Euclidean rotations or taste transformations, but only under the 
lattice symmetry where taste and Euclidean rotations by 90° are performed simultaneously. 

In Refs. and the leading order (LO) chiral Lagrangian for light-light pseu- 
doscalars is found for the one-flavor and many-flavor cases, respectively. There the only 
available Lorentz 4- vector in the chiral theory is the partial derivative operator 5^. To con- 
struct a chiral operator with the symmetries of operators requires at least two derivatives 
and two explicit taste indices. Such operators are O(p^a^), where the factors of lattice spac- 
ing a come from the explicit violation of taste symmetry, and the factors of p come from the 
derivatives. These operators then only occur at next-to-leading-order (NLO) in the chiral 
Lagrangian, unlike the representatives of which are LO. The representatives are 
therefore omitted in Refs. jlll. If^j]. although they are included in Ref [2^, which works to 
NLO. 

In the heavy-light case, however, there are two 4-vectors available in the chiral theory 
whose presence does not raise the order of the operator: the heavy quark 4-velocity 
and the light quark gamma matrix 7^. Therefore the chiral representatives of the 
operators contribute at the same order as the representatives of the C^. For heavy-light 
decay constants this is NLO — such operators enter into analytic terms in the calculation 
presented here, but not in one-loop diagrams. 
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III. THE SXPT LAGRANGIAN WITH HEAVY-LIGHT MESONS 



A. Leading order continuum theory 



We first discuss standard XPT at LO for lieavy-liglit mesons in tfie continuum 3lLl32l.l33, 
Is^ . More precisely in tfie current development, tfiis is tfie LO cliiral representative of the 
dimension 4 part of the Symanzik Lagrangian, Cf^, Eq. Due to the heavy quark spin 
symmetry in the static hmit, the heavy vector and pseudoscalar mesons are incorporated 
into the following field, which destroys a heavy-light meson 

Ha = ^ b'B;^ + ^1,B^] , (8) 

where v is the meson's velocity, and a is the combined light quark flavor-taste index. The 
conjugate field creates a heavy-light meson 

Ha ^ loHho = [YBl: + ^75^1] ^ , (9) 

We use B to denote a generic pseudoscalar heavy meson and B* to denote the corresponding 
vector meson, but note that the current formalism will also apply to D and D* mesons, 
although at a decreased level of accuracy in the l/mg expansion. 

Under the SU{2) heavy quark spin symmetry, the heavy-light field transforms as 

H ^ SH , 

H HS^ , (10) 
with S G SU{2), while under the SU{4:n)L x SU^An)^ chiral symmetry, 

H HV^ , 

H ^ VH , (11) 

with U G SU{An), and we are keeping the flavor-taste index implicit. 

We use the term "pion" generically to denote any of the light pseudoscalar mesons. The 
field describing the pions is S = exp [«$//], where / is the tree-level pion decay constant 
(/ ~ = 131 MeV). In terms involving the heavy-lights, we also need a = = 
exp[i$/2/]. These fields are singlets under the heavy-quark spin symmetry, while under 
5'f/(4n)i X SU{An)ji we have 

S ^ LJ:R^ , ^ i^S^L^ , (12) 

a LaV^ = UaR^ , Ra^U^ = Va^L^ , (13) 



with L G SU{An)L and R E SU^in)^. Equation (fT^ defines U, which is a function of the 
pion fields and the coordinates, as well as of the global transformations L and R. Chiral 
operators are formed from combinations of H, H, a, and derivatives, as well other matrices 
such as the light quark mass matrix or taste matrices, so that they transform in the same 
manner under the combined symmetry group as the underlying operators in the Symanzik 
action. As usual, this is accomplished by promoting any symmetry-violating factors in the 
Symanzik Lagrangian to spurious and choosing the transformation properties of the spurious 
to make the Lagrangian invariant. 

For n KS flavors, S = exp(i$//) is a 4n x 4n matrix, and $ is given by: 



u 


71+ 


K+ 




D 




K~ 




S 



\ 



V 



(14) 



where U = X]==i U-eT-e-, etc., with the Hermitian taste generators given by 



(15) 



As in Ref. |l3| the taste matrices are with = ^^^5. The 4x4 identity matrix is 
^/ = I. In this paper, we define ^^i, = (1/2) [^^, ^,^] (rather than simply to emphasize 
that terms with ^ = v are never included; in Eq. p5|l we take fi < u only. 

The component fields of the diagonal (flavor-neutral) elements (f/n, D^, etc.) are real; 
the other (charged) fields are complex (vri, K~, etc.), such that $ is Hermitian. Here the 
n = 3 portion of $ is shown explicitly. The mass matrix is the An x 4n matrix 



/ 



M 



mj 
mJ 
mJ 



\ 



\ 



(16) 



where again, the portion shown is for the n = 3 case. 

It is convenient to define objects involving the a field that transform only with U and 
The two possibilities with a single derivative are 



(17) 
(18) 
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transforms like a vector field under the SU{4n)L x SU{4n)ji cliiral symmetry and, when 
combined with the derivative, can form a covariant derivative acting on the heavy-light field 
or its conjugate: 

0;H), = DfH, ^ d,Ha - iVfHb , (19) 

with implicit sums over repeated indices. The covariant derivatives and transform under 
the chiral symmetry as 

A^ ^ UA^U^ . (20) 
We then write the (Minkowski space) continuum LO chiral Lagrangian as 

'^LO,cont = '^LO^'cont + ^1 ' (21) 

where -^Lo'cont standard LO light meson Lagrangian in the continuum, and Ci is the 

leading term involving the heavy-lights. As explained in Sec. 1111 Bl lattice corrections affect 
the heavy-light terms only at NLO, so there is no need to specify "cont" on Ci. We have 

^IZont = Y Tr(5,S9-St) + Tr(A^S + ^St) - J^^Uj + Dj + Sj + . . .f (22) 



and ^ ^ 

£i = -t Tt(HHvD) + Tr(M7'^75A^) . (23) 

Tr means the complete trace over flavor-taste indices and, where relevant, Dirac indices. 
Since H and H always appear together in the Lagrangian, we treat HH as a matrix in 
fiavor-taste space: {HH)ab = HaHf,. Here and below, the covariant derivative D acts only 
on the field (in this case H) immediately preceding it; similarly D acts only on the field 
immediately following it. As in Ref. it is helpful for analyzing the quark flow to leave 
the anomaly {ml) term explicit in Eq. (|22j) and work with the diagonal fields U, D, S, . . . , 
rather than the physical ones ttq, tl- ... At the end of the calculation, we can take rriQ oo 
and return to the physical basis |35l | . 

Although fundamentally we are interested in the chiral effective theory for a Euclidean 
lattice theory, we find it more convenient to write the Lagrangian in Minkowski space to 

11 



make direct contact with the extensive continuum hterature. A Wick rotation to define the 
Euchdean Green's functions is imphcit everywhere. 

For the heavy-hght decay constants, we also need the chiral representative of the axial 
heavy-light current. Alternatively, one can work with the left-handed current, which has the 
advantage that it transforms more simply under chiral transformations and in addition gives 
corresponding information about the vector current, useful for semileptonic form factors. 
The left-handed current that destroys a heavy-light meson of flavor-taste b is j^'^, which at 
LO takes the form 3^ 

to = ^M^'i^-^^)Hy^^'^ (24) 
where A*-*-* is a constant vector that flxes the flavor-taste: {X^''^)a = Sab, and tr^ is a trace on 
Dirac indices only. In QCD, the decay constant fs^ is deflned by the matrix element 

{0\r'\Baiv))=tfB^mBX5ab, (25) 

where relativistic normalization of the state \Ba{v)) is assumed. At LO in the heavy-light 
chiral theory, j'^q = inv^Bf,, which gives f}^^ = k / ■ Recall that the factor ^JtubI 
arises from the differences in normalizations between relativistic and non-relativistic states. 



B. Power counting 

Before considering discretization errors and higher order corrections to Eqs. (j2H) and (j24j) . 
we discuss the power counting assumed in this paper. The staggered chiral Lagrangian with 
only light mesons is a joint expansion in the light quark mass and the lattice spacing a^. 
In that case, and are taken to be of the same order in the expansion jlll. u^. Since 
oc vnq and the momentum of external light mesons, p^, is assumed to be of order m^, we 
have ~ ~ ~ in the power counting. 

Upon including B mesons, we have an additional expansions in the inverse of the heavy 
quark mass mg, and in the B meson's residual momentum, fc, which we take to be of the 
same order as We work to leading order in l/mg only. 

The continuum LO chiral Lagrangian, Eq. ()2H). is therefore Oik ~ y/fng) in the heavy- 
meson flelds and 0{mq, a^) in the light meson flelds. In each case these are the lowest 
order terms allowed by the continuum symmetries.'^ One-loop diagrams with vertices from 

^ The 0{k'^) heavy quark mass term is absorbed by measm'ing heavy-hght energies relative to the mass of 
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Eqs. ()2ip and ()24|) and propagators from Eq. ()21|) will give O(m^) corrections to physical 
quantities such as fs^ or semileptonic form factors. Since taste- violating terms of O(a^) in 
the pion Lagrangian of SXPT are the same order as those already included in Eq. (j221), we 
will need to include such terms in order to calculate the one-loop corrections. 

Terms that are higher order than those in Eqs. (j21|) and (j24p are irrelevant in the one- 
loop diagrams. They can, however, make analytic (tree-level) contributions at this (or, in 
principle, lower) order. In particular, for fs^ we will need to consider 0{^mq) and 0{mq) 
corrections to the existing heavy-light terms in Eqs. ()2H) and (j211). Such terms in the heavy- 
light Lagrangian, Eq. are of order niq ~ a? and kmq ~ ka^; those in the current, Eq. ()24j) . 
are of order ~ o?. These corrections are generated by inserting, into heavy-light chiral 
operators, spurious representing the light-quark mass term in Cf [Eq. ©] and the 0{a^) 
taste- violating terms, [Eq. Q]. There are also invariant higher order operators that are 
present even in the continuum and chiral limits. These are generated by inserting additional 



derivative operators, i.e., additional powers of k. As explained in Sec. IIII D[ such terms do 
not affect the decay constants to the order we are working, and we do not attempt to catalog 
them completely. We also do not need to consider here higher order [0{mg ~ m^a^ ~ a^)] 
corrections to the pion Lagrangian. Such corrections cannot give analytic contributions to 
heavy-light decay constants. They are, however, relevant to purely light-light physics, as well 
as to semileptonic form factors of heav y-lig ht mesons, which involve a light-light particle. 
All such terms have been found in Ref. (8fll |. 

Summarizing the discussion of this subsection, we write our complete chiral Lagrangian 

as 

C = Ci^o + 'C^Lo , (26) 

^ ^pions . ^ ^pions , ^pions , ^ ^07^ 

f-LO - ^LO + - ''-LO.cont + ^LO,a:2 + , \^<) 

-^^NLO = + (28) 
£2 = (^2,k + £2,™ + £^a2 + ^2,0^ 1 (29) 
£3 = ^^■i^k + £3,™ + -^^a^ + -^^02 • (30) 

Here, £1, £21 and £3 involve the heavy-light fields and are of order fc, k"^ ~ ~ a^, 
and k^ ~ kmq ~ ka^ , respectively. The subscripts /c, m, or 0? indicate terms that involve 



the heavy meson. 
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derivatives only, mass spurions, or taste-violating spurions. The superscripts A and B stand 
for the type of taste- violating terms, Eqs. (jH) and (0). 
Similarly we have for the left-handed current: 

J'^' = J^i + J^to, (31) 

J^to = J^'' + Jt (32) 

Ji"' = Jtf (33) 

J^' = Jtk + & + ji^A + li^B ■ (34) 

Note that insertions of the mass and a? spurions are always at least two orders higher in 
kr^y/fn^, so do not appear in jf'''. 

The continuum LO terms Ci, /^Lo'^cont) Jlo given in Eqs. (j^ . (^Hj) and (j^ . The 
remaining contributions, Lagrangian terms 'C^^^'^, 1^2-, and £3, and current terms jf'^ and 
22^^ are discussed in the following subsections. 

C. LO discretization effects 

The LO finite lattice spacing correction, C^^^p.-, to the pion Lagrangian is known from 



previous work 



11 



3 



We have 



Vs = CiTr(ef Ee^")st) + C3^$^[Tr(e)se(")S) + /^.c.] 

V 

V V 

+ ^5y^ 5^ Tr(e^S) Tr(e^St) + Csa^ Tr(d')s) Tr(^(?St) . (35) 

The explicit 4n x 4n matrices in Eq. ()35p are defined by 

(^^"^).,=e.5.. > (36) 

with i and j the SU{n) (light quark) flavor indices, and the 4x4 taste matrix, as in 
Eq. (fT3j). The matrices and ^^g'' are defined similarly. 

There are no corrections at LO to Ci, Eq. (|^. or to the left-handed current, Eq. 
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D. Invariant NLO corrections; T symmetry 

Even in the chiral and continuum limit, there are invariant terms that can be constructed 
by adding more derivatives to the lowest-order heavy-light Lagrangian £1, or current j^Q. 
To our knowledge, a complete catalog of such terms through 0{k^) does not appear in the 
continuum literature. These higher derivative terms are in fact irrelevant for the heavy-light 
decay constant to the order we are working: Derivatives acting on the light-light fields E 
or a do not contribute to tree diagrams because there are no external pions in the matrix 
element; while terms with a derivative acting on the heavy-light field in the current, or 
more than one derivative acting on a heavy-light field in the Lagrangian, give a vanishing 
contribution on shell {k — 0).^ 

For semileptonic form factors, e.g., those for B ^ ti, higher derivative terms do produce 
analytic corrections at this order. However, the functional form that such contributions 
can take is rather limited: In the chiral limit, one can only get corrections to form factors 
proportional to v-p or {v-py, where p is the pion 4-momentum. Away from the chiral limit, 
there will be additional corrections proportional to p^, but these merely duplicate, on-shell, 
the effects of mass or spurious. 

We thus simply list here a few representative operators that appear in C2,k, C,z,k-, iff, 
and J2^f : 

C2,k = 5^ Tr {{v ■ DHH - HHv ■ D) -f^-f.^) + ^ ^r (h{v ■ DfH) +... (37) 

X X 

C^,k = || 1^ {hH^M''^ ■ Df^''^ + |i Tr (hH^^, vDvAj+... (38) 

X X 

Jii - ^tr,(7''(l-75)^-Di/)atA('') + ^tr,(7'^(l-75)i/)^-AatA(^) + ... (39) 

X X 
X 

M^"" {^-l^)H)v ■ Dv • AatA(^) + . . . (40) 

X 

where the constants e^, Kj are taken to be real and dimensionless, A^ is the chiral scale, and 

D,A^ = d,A^-i[Y,,A^] . (41) 

Terms with a single derivative in the Lagrangian do contribute on shell through wave-function renormal- 
ization. 
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The use of time reversal (T) symmetry and the requirement of Hermiticity for the La- 
grangian can be used to ehminate several other candidate operators from Eqs. (jHTf) through 
(|4Up ■ Because we will also need to use the consequences of T invariance extensively in 
Sees. 1111 El to llV Bl we briefly review how it acts on relevant objects, following Ref. 0. Let 
p be the parity reflection of a 4- vector p; in other words p^ = Pf_i- Then time reversal acts 
as follows 



X —X 
V ^ V 



^{x) - 


-$(-a; 


) 


- 






a{x) - 


a{-x) 




A'^(x) - 


A^(-x) 




V^(x) - 






Df" - 






H^ix) - 


TH\- 


x)T'^ 


H\x) - 


th\- 


-x)T^^ 



(42) 

where the Dirac matrix T has the property T'-y^T~^ = 7*, we show the v dependence of H 
explicitly, and we have used the anti-unitary nature of T. 

Requiring T invariance, we can for example eliminate contributions to the Lagrangian 
such as 

Tr (m7^75^; ■ :D A^) . (43) 

Add,.g a face, of . to .aUe this T-even. as was do.e ,n Ref. Q . does not faet save tMs 
term because it is then anti-Hermitian. 



E. Quark mass corrections at NLO 

Here we discuss operators induced by single insertions of the light quark mass spurious in 
£f^. These produce corrections to the Lagrangian, namely £2,™ and £3,^, and corrections 
to the current, J2 m- course, such terms are also present at NLO in the continuum, and 
the result is more or less standard 0, IsTj. However, since many of the contributions to 
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the discretization corrections {^,202^^20^^ '^sa^' '^fa^' ^2aP-A^ ^"^^ ^2 a2 b) follow the same 

carefully 



36 



37| 



pattern (see Sees. llVXl and UVBjl . it is worthwhile to examine the mass corrections__ 
first. In addition, we find some modifications to the terms written down in Refs. 

We start by writing down combinations of the light quark mass spurious that transform 
only with U or These are 

A1± = 1 [aMa ± a^Ma^) , (44) 

which transform as — > UA^^'^U^ under the chiral symmetry. Under parity, M."^ — *• 
iAl''^, since $ — $ and cr cr''"; while both M."^ and }A~ are even under time reversal. 
Because tTo{HH'~f^) = 0, the only parity and chiral invariants we can construct with no 
derivatives are TY{HHAi~^) and Tt{HH) Tr(A^^). We thus have 

£2,m = 2Ai Tr {HHM+) + 2X[ Tr (HH) Tr {M+) . (45) 



For Cs^m, we start with the lists of terms given in Eq. (16) of Ref. |3a] and Eq. (7) of 
Ref. j37|. Note, however, that a term like Tr (^HH'j^'j^ A^Ai^) is not Hermitian, because 
the Hermitian conjugate will flip the order of and A'^. To make it Hermitian we can 
write it in one of two ways 

zTr (M7^75[A^Ai+]) . 

The second term is not invariant under time-reversal and can be dropped. For corre- 
sponding terms with A4~ , it is the anticommutator that requires a factor of i, since 
is anti-Hermitian. Then time- reversal invariance requires that we build the traces 
using the commutator [A^,Ai'] only. Terms with separately traced operators, such as 
Tr (if if 7^75) Tr(A'^Ai+), do not need to be rewritten with commutators or anticommuta- 
tors since they are Hermitian already. 

Terms involving a derivative acting on the heavy-light fields also require some thought. 
Hermitian combinations are 

i{d;hh-hhd^ 

D;HH + HHD^ . 

However, we will not need to include the second combination, which is a total derivative. 
When combined a factor of A^^, integration by parts can be used to put the derivative on 
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the 0" or a'^ fields. Such terms can then be rewritten in terms of A, using iD^a"^ = —A^^cr^ 



or iD^a = A'^cr 



36|. 



Putting together the above discussion, we have 

+ik2 Tr (hHvD-vDHH^ Tr(A^+) 

+ A;3 Tr {HH^^^,{A^^, M+}) + k, Tr (HH^^^.A^) Tt{M+) 

+ k^ Tr (HH^^j^) Tr (A^A^+) + k^ Tr (m7/.[A^ M']) . (46) 

Note that H and H must be next to each other in these traces by heavy quark spin symmetry. 
The term that involves A4~ will not contribute to the decay constants or to the form factors 
for semileptonic decays, but we include it here to complete the list of mass-dependent terms. 
For the current, we have: 

+p^tT^{Y{l-l5)H)M-a^X^'^ +p^tTo{Y{l-l5)H)cx^X'^''>TT{M~) . (47) 

It is not hard to argue that the heavy quark spin symmetry forces the factor 7^(1 — 75) to 
appear next to the field H; it could not, for example, be replaced with v^{l — 75). To show 
this, one can replace 7^(1 — 75) in the QCD current by a spurion, and give it appropriate 
transformation properties under heavy quark spin symmetry and Lorentz transformations 
to make the current invariant. A similar argument will apply even for the more complicated 
terms in the current described below. 



IV. TASTE BREAKING WITH HEAVY-LIGHT MESONS 

The final needed contributions to the chiral Lagrangian and current are 0{a'^) terms 
involving heavy-light fields: £^^2, £|„2, £^„2, £^^2 [Eqs. and and j2'a2,A' ^v^s 

[Eq. ()34|) ]. These are constructed by starting with the dimension-6 terms in the Symanzik 
Lagrangian, and [Eqs. (jlj) and (0)]. We promote the two explicit taste matrices in 
each of these operators to spurion fields, which are assigned transformation laws such that 
the operators are invariant under the chiral SU{An)L x SU{An)R symmetry. 

The heavy quark spin symmetry requires that we only have the combination HH in the 
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Lagrangian, and 7'^(1 — 75)if in the current. Under the chiral symmetry, these transform as 

HH VHHV^ 

When doing the spurion analysis, we therefore need to combine the taste matrices (treated 
as spurion fields), a, and cr^ into operators Oi that obey: 

Oi UOiV^. (48) 

We then construct chiral (and heavy-quark spin) invariants using these spurious and the 
heavy-light fields. At 0{a'^), exactly two taste spurious will appear in each term in the 
Lagrangian and the current. In C^^2 and >C^^2, which are 0{ka?), a derivative term will 
also be needed. Either or can be used. 

A large class of contributions to £-^^2 , C.'^a? ' ^'^'^ 3 -i t? a "^^^ obtained by first combining 
both taste spurious into a single operator that transforms like Eq. (j48|) . As described in 
Sec. IIV Al all such operators may be easily found from known results using a trick. These 
operators are the same order in our chiral power counting and have the same transformation 
properties as the quark mass terms, }A^. Thus, we can simply replace the M."^ in Eqs. 
fliH|) and (jTTj) with such two-taste-spurion operators to find all terms in >C^^2, i^^a?i and 
^2 0? A which the two spurious combine together into a single operator. It is then not 
hard to catalog the remaining contributions, which are ones in which the two spurious are 
interspersed between HH and A in a single trace, or are separately traced with HH or A. 

Determining the chiral representatives of works similarly. However each chiral invari- 
ant constructed from the spurious will generate several terms in £^^2, -^^02, and due 
to various ways of introducing the appropriate Lorentz indices using w^, 7^, D^, or A^. 



A. Discretization errors at NLO: Operators from 

To find combinations of two taste spurious from that transform like Eq. pH|l . we can 
start with the invariant light-meson operators that result from the spurion analysis. These 
are the eight operators given in Eq. dSSl) QQ. From these operators, we can generate all 
ones that transform like Eq. ()48|) by simply replacing a factor of S or with, respectively, 
(T^ or (o"^)^, permuting the trace so there is a factor of a or at both beginning and end of 
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the term, and removing the trace. For example, 

Tr(d"^SteJ")s) ^ a^-^St^J-V = a^-^tat^^-V , (49) 

where the final form shows explicitly how the result is composed of the produce of two taste 
spurious, each transforming like Eq. ()48|) . We call this process "opening up" a S (or S^). 
It automatically creates an operator that has the desired transformation property when the 
spurious transform as required. Each operator in Eq. thereby generates two operators, 
depending on whether one opens up a S or a As was done for the mass term in Eq. 
the two operators can be joined to make Hermitian or anti-Hermitian combinations, or, 
equivalently, even or odd parity combinations. This gives us the following operators that 
transform under chiral rotations by Eq. ()48p : 





= (^d"^Ste("V±p.c.) 




























V 


Kif^)Tr(d?S)±p.c.^ 


^8 


V 


K£^^)Tr(d?St)±p.a 



where p.c. stands for parity conjugate: for example cjp.c. = o"T. We have followed the number- 
ing scheme in Eq. (j35|l . with the exceptions that the operators multiplied by C2V and C^y 
in Eq. become C^'^ and C^'^, while those multiplied by C2A and C^a become O^'^ 
and O^'^ . We have also not bothered to keep the conventional factors of 1/2 and 1/4 from 
Eq. fl35|) in our definitions here. Operators in Eq. ()50p with no traces (coming from ones 
in Eq. (jH^j) with one trace) correspond to multiplying two taste spurious transforming like 
Eq. ()48p: operators with a trace (coming from ones in Eq. ()35p with two traces) correspond 
to tracing one taste spurion and then multiplying by the other. 
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Making the replacement A^^ — > O^'^ in Eq. ()45p . we then get 

^ta^ = Tr {HHOt'^) + Ki, Tr (M) Tr(0,^'+) I (51) 



fc=i 

C^^2 is more comphcated, since we can use one factor of A'^. Because transforms 
hke HH , there is now an additional possibility of interspersing single spurious between A'^ 
and HH. A complet e sp urion analysis is therefore required for such terms. We can follow 
Appendix B of Ref. ll|, with the simple modifications 13] that we keep only "odd-odd" 
operators and make the spurious diagonal in flavor, as in Eq. (jHEI)- We then prepend and 
append appropriate factors of a and o"^ to each spurion to make it into an operator that 
transforms like Eq. PHj) . The result of this analysis is that a given chiral operator of this 
kind in C'^^p. will involve a pair of spurion operators, which we call and (/c = 1, . . . , 8): 

















— (Qi )p.c. 








[VxP], 










^Qi 




[S X 




[Px V], 


[T X -> 


Qf 


^(n) 


0^ 


^Qi 




[S X 




[Px A], 


[TxA] 


Qt 




Of 


^Qi 








[5x\/], 


[PxV]^ 


Qi 


^(n) 




= {Q5 )p-c. 








[\/xT], 


[AxT]^ 


Qi- 


■ An) t 


Qi 


= (Qe )p.c. 








[KxT], 


[AxT]^ 


Qf- 


• An) t 


Qi 


^Qi 








[S^A], [P X A] ^ 


Qi 




Qi 


= (^8 )p-c. 





where we show the terms in [Eq. (^] that generate the given spurion operators. The 
terms [V x S] and [A x S] in the Symanzik action do not appear in this list since the 
corresponding spurion operators, such as a^^V^, are trivial. We have tried to make the 
numbering system in Eq. ()52p correspond to that in Eqs. (jHSj) and ()50p as much as possible. 
Note that = Q2 but Qf ^ Q^, and similarly for = Qf, Qf = Qf, and = Qi- 
Operators obtained by parity conjugating both and are not included since they will 
appear automatically when we demand parity invariance of the chiral Lagrangian. 

We can now construct ^33^2 by first replacing Ai"^ — C^'^ in Eq. (jlEI), and then adding 
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on new operators that intersperse and between HH and A^. The result is: 

k=l I 

+ic^fc Tr (hHvD-vDHH^ Tr(C^'+) 

+4, Tr (m7,75{A^ Of +}) + c^, Tr {HH^,^,A^) Tr(0,^'+) 
+c^, Tr (m7M75) Tr(A^O,^'+) + c^, Tr (hH^,[A>^, Otn) 
+c^,(Tr (m7^75gf A'^Q^) +p.c.) 

+4, ( Tr {HH^,^,Q^) Tr (A'^Q^) + p.c.) | 

W Yl 4k{T^T{HH^,QtA^Qt)+p.c.) 

fc=2, 5,7,8 

E ( (HH^.Qi) Tr (A'^Q^) + p.c.) . (53) 

fc=l,2,6,7 

Here, any taste indices in Q'^ and should be contracted; for example one should sum 
over A and i/ in a term involving Q^, Qq.^ Note that only a subset of the are allowed in 
the terms with coefficients Cg^ and c^ofc' which involve 7^ as opposed to 7/^75. The parity 
conjugation introduces a minus sign, so the terms can only be Hermitian if = {cgf^ 
term) or = (c^ofc term). The forbidden terms are anti-Hermitian; adding a factor 
of i to make them Hermitian would violate T invariance. 

The reader may wonder whether additional operators in C^^2 might be constructed us- 
ing various Qf and Qf in place of the Cf ^. The answer is no, since the correspond- 
ing terms would contain either the product QfQf or QfTi^Qf) (along with appropriate 
parity-conjugated terms). In these circumstances, it is easy to check from Eq. (j^^ that one 
reproduces the Cf ^. 

The converse question is why the Qf, Qf cannot be found by simply opening up a 
second S or S''" in the Of*. The reason is that certain allowed combinations of spurious — 
namely Q21Q2 Qii Qi — have been eliminated from Vs and hence from Cf ^ because 
their products are trivial and their traces vanish. Such combinations can only be found by 
returning to the original spurion analysis. 

^ Contributions with \ — v are automatically omitted since we use the definition ^aj/ — (1/2)[Ca, ■Ci']- 
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)p.c. 






^(n) 






)p.c. 






■ A^) 






)p.c. 






• t 




= QiiA ' 







For the current at this order we cannot use A^, so the replacement — > 0^'^ in 
Eq. (P7|) is all that is needed. We have 

^i,A = Ej^i tr, (7^(1 - 75)^^) Of VtA(^) + 4, tr, (7'^(1 - 75)^^) ^tA^") Tr(0,^'+) 
fc=i I 

+^3^. tr, (7'^(1 - Ot'-a^X^'^ + r^, tr, (7'^(1 - 75)^^) aU^^) Tr(0,^'-) | . (54) 

B. Discretization errors at NLO: Operators from 

The terms which arise from breaking rotation symmetry at this order can be determined 
through a spurion analysis that follows that of Sec. A3 in Ref. js^. The following pairs of 
spurious can appear: 

[V,xT,], [A, XT,] 

[V, X T,] , [A, X T,] ^ Q^^, = ^a^ilVt , Q^, ^ Q^, , (55) 

where we show the terms in [Eq. (jSj)] that generate the given spurion operators. The 
index /i is singled out in the names Q^j because it will be repeated 4 times in terms in C^^2, 
thereby breaking separate continuous rotations and taste symmetries, and leaving only the 
lattice symmetry of joint 90° rotations. We do not include the index A in the names Q^i 
and because A will appear only twice, and a sum over A will be implied in any terms 
with Q^^i, Q^^i or Q^4, Qj^^4- The spurious Q^j (j = 1,2,3,4) are in fact the same as Qq, 
Q^, Qg and Qf, respectively. 

We can also construct operators, analogous to the Of'^, out of either Q^jQ^j ± p.c. or 
Qjj.Tr(Qj^.) ±P-C- These are 

O^^t = (aei"V)Tr(ei'^)st)±p.c. , 

O^t = Ki?a)Tr(e£)st)±p.c.. (56) 

The numbering here corresponds to that in Eq. Q^4, Qj^4 produce no operators of this 
kind. Note that the repeated index n is not summed over in Eq. aside from this these 
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operators are the same as operators C^'^, Of'"^, and C^'^, respectively. Like the the 
O^'f^ transform according to Eq. (gS)). 

To obtain terms that correspond to C^, we can multiply either the O^'^ or the Q^^, 
Q^,k by two additional four- vectors with index /i, and sum over /i. For >Cf^2, we can employ 
Vfj_ or 7^. We find: 

= E Tr(MOf '+) + K.^.t;,^;^ Tr(M) Tr(0^,+) 1 (57) 

We have used that fact Ti{HH^^ oc T^i{HH)v^ to eliminate other possibilities. Terms with 
Q^fc, Q^fc are redundant and can be rewritten in terms of O^'^ . 

For C^^2i a single factor of Z^^ or is also allowed. The enumeration then gets quite 
complicated, because Q^^, Q^^k can appear in the terms, and, in addition, various familiar 
simplifications cannot be used. For example, manipulations that follow from Hi/) = —H will 
not apply when the summation index ^ m. = also appears in the taste-violating 

operators. To make the equations more manageable, we write, 

Cl^. = + Cl;^ . (58) 
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We then obtain: 

fj, k=l I 

+cj, IV (m7m75{A^ + cj, Tr (HH^,^,A^) IV(Ojf ) 

+cg, Tr (HHj.j,) Tr(A'^Oj)+) + cj, Tr (hHj,[A^, OJ-]) 
+ic^, TV (hHvDO^J - vDHH 

+^C3^,^;XTr (hHvD-vDHh) Tr(Of'+) 

+cj, TV (S^i^7.75{A^ + cfo,, ^^^'^ TV (HH^.^.A'^) TV(Of 

+cfi,, r;,^'^ Tr (M7.75) ^(A'^O^;^) + cf^,, Tr (:ffi^7.[A^ ^J^]) 
+cf3,, Tr (:ffi/7^75{^-A, + cf^,, Tr (:ffi/7M75 ^-A) Tr(0^^)+) 

+cf5,, TV (ra7M75) TrivA O^J) + cfe,,^;'^ TV (hH^,^^,{A\ '+}) 
+cf7,,^'^ Tr (:ffi/7..75) Tr [A^Of^) + cf^,,^'^ Tr {HHj,^j,A-^) Tr (o^f ) 

+4,fe^" Tr (m7^.{A^ Oj'^}) + Tr (HH^,^) Tr (A'^Of ),-) | (59) 

A sum over the repeated index u is imphed, and 7^,^ = (1/2) [7^, 7^,]. We have used the fact 
that Tr(C^^~) vanishes for aU k. 
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For C^'^, we find: 



B,Q 
3,a2 



. fc=i >- 



+4,fc(Tr {HH^,^,Ql,) Tr (A^'QJ,) +p.c 

+4,,t;X(Tr {HH^.^.Ql,) Tr (A^QJ,) +p.c. 
+4,,i;'^(Tr (M7M75QS.t;-Ag2,) +p.c.) 
+4,,t;'^(Tr [HH^,^,Ql,) Tr (t;-AgJ,) +p.c. 
+4,,T;'^f Tr (M7M.75g?,,A^g2,) +p.c. 

Tr {HH^,,^,Ql,) Tr (A^g?,.) +p.c. 
+ E 4,fc ( Tr (M7MQS.A^gJ,,) + p.c) 

E 4,. ( Tr (HH^.Ql,) Tr (A^gJ,,) + p.c 



/c=2,3,4 



k=lA 



Tr(m7.QM)Tr(A^gSfc)+p.c. 
+43,it;'^( Tr (M7m.Q?,i^^^?,i) + P-^ 

+ E 

fc=2,3 



4,,^;^ Tr {HH^,,Ql,) Tr (A^gJ^,) +p.c 



(60) 



The additional implicit taste index A that appears in g^i, Q^i or g^4, g^4 should be 
summed. As in Eq. ()53p. only a subset of the Q^i^ are allowed by Hermiticity and T 
invariance in terms where the parity conjugate introduces a minus sign. Since A is an axial 
vector, these are the terms without an explicit 75. 

Note that we have raised one space-time /i index in each term in Eqs. (j59p and 
because we are using Minkowski space conventions. We also emphasize that the 7^ here, 
and elsewhere in the chiral Lagrangian, are always Minkowski ones. The taste /i index on 



O^'t and gf/fc, gf/fc is not affected by going to Minkowski space, and taste matrices ^„ obey 



' fj,,ki ^fj.,k 



Euclidean conventions. 



26 



Proceeding to the current, we have: 
32:a^,B = Mtr,(7^(l -75)i^)Of.VA('') + rj,tr,(7'^(l -75)i^)atA 

k=l I 

+ri,, tr, il'il - 1,)H) Ojr^^AW + tr, {^{1 " l^W a^X^'^ Tr(Of 
+ E (^5"!. tr, (T'^ll - 75)^^) t;.t;^a^,VtA(^) 

+ri^, tr, (7'^(1 - 75)i/) ^U^^.t;'^ MO^) 
+r,^, tr, (7'^(1 - l,)H) v,v^ 0% a^X^'''^ 
+ri^, tio iril - l,)H) a^\^'\,v^ Tr(0,^r) 

tr, (7'^(1 - 75)i^7^) ^^.aY^^^^'^ 
+rfo.. tr, (7^(1 - 75)^7^) <y^\^'^v, Tr(0^;,+) 
+rfi,. tr, (7^(1 - 75)i^7^) v,0^^;,a^\^'^ 

+<.tr, (7^(1- 75)^7^) ^U(%.Tr(0;;-)) I . (61) 

Here the repeated index /i is not summed over. 

Despite the complexity of Eqs. through (jHH), the effect of these terms on the heavy- 
hght decay constant is extremely simple. At this order, the terms can contribute only at 
tree level, so we can set cr = 1 and = 0. Each term then either vanishes outright, is 
nonvanishing but does not contribute to /s, or reduces to an additive constant in fs that 
is proportional to a?. In particular, there are no rotation- violating contributions to until 
higher order. All these 0? analytic terms thus combine and lead merely to the presence of 
one additional unknown parameter in SXPT chiral fits, compared to continuum- like chiral 
fits. 



V. CHIRAL LOGARITHMS IN Jb AT ONE LOOP 

Before calculating the one-loop diagrams, we need to write down the propagators of the 
heavy-light and pion fields, as well as the vertices coupling them. We use the notation 
{V'xK^) to denote the Minkowski space propagator of fields ip and x with momentum k. 
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Then, from the LO heavy-hght Lagrangian £i, Eq. ()23|). we have 



B;.B:l}(t) = -"^1^^;;^"^ . (63) 

Here a, b indicate the flavor-taste of the hght quarks. Lower case Latin indices from the 
beginning of the alphabet {a,b,c,d . . .) will always serve as such flavor-taste indices; in the 
full QCD case they run over 4Nf values, where Nf is the number of light sea-quark flavors. 

Since taste violations do not appear in Ci, Eqs. (lU^ and (jU^ are the same as in the 
continuum theory js^ , except that flavor-taste indices replace pure flavor indices. Similarly, 
the BB*7T vertex looks identical to that in Ref. [s^ aside from the redefinition of the indices 
and a factor of 1/2 in the present case coming from our normalization of the taste generators, 
Eq. ()15|1 . The term in the interaction Lagrangian that gives this vertex is: 

y Wai?. -i?Ii?;.)5'^<^'^.a, (64) 

where repeated indices are summed. The other needed vertex comes from the Bnir term in 
the expansion of the LO current, Eq. ()24j) . We have: 

ti = tnv''(^B,-^B,<i>A + ---^ , (65) 

where repeated indices are again summed and ■ ■ ■ represents terms involving B* or other 
numbers of pions, which do not contribute to the decay constant at this order. 

When necessary, each flavor-taste index can be replaced by a pair of indices representing 
flavor and taste separately. We use Latin indices in the middle of the alphabet {i, j, ■ ■ ■) as 
pure flavor indices, which take on the values 1,2,..., Nf in full QCD. Greek indices at the 
beginning of the alphabet (a, /5, 7, . . . ) will be used for quark taste indices, running from 1 
to 4. Thus we can replace a ia and write, for example. 

For light mesons, the LO Lagrangian C^^^, Eq. includes both the continuum terms 
'^Lo'cont leading corrections As in Ref. Q], we treat explicitly the "hair- 

pin" terms arising from the mg contribution to Eq. and the C2V, C2A, C^v, and C^a 
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contributions to Eq. by separating out the disconnected parts of meson propagators. 
The remainder of C^q'^'' determines the connected pion propagators: 

I. J conn p — ij = ~r 

where S is one of the 16 meson tastes [as defined after Eq. ()14|) ]. and rriij^s is the tree-level 
mass of a taste-S meson composed of quarks of flavor i and j: 

= = ii{mi + TUj) + a^Ac. (68) 

Here Ac is the taste splitting, which can be expressed in terms of Ci, C3, C4 and Cg in 
Eq. dSSD [13]. 

Because of the residual 5*0(4) taste symmetry [ll| at this order, the mesons within a 
given taste multiplet (P, V, T, A, or /) are degenerate in mass. When it is not important 
to specify the particular member of the multiplet in question, we will usually just name the 
multiplet, for example: 

fn^jy = /^("^i + "^j) + ct^Ay. (69) 

Since the heavy-light propagators are most simply written with flavor-taste indices, as in 
Eqs. (jU^ and (jU^ . it is convenient for current purposes to rewrite Eq. (p7j) in fiavor-taste 
notation also: 

{*,Av} (P) - (P) = E T'r?^t • 

J conn I. J conn •^l^ p — iTi'^n \ i^t 

where are the 16 taste generators, Eq. (fT^ . 

For flavor- charged pions (z 7^ j), the complete propagators are just the connected propa- 
gators in Eq. ()67p or (f7n|) . However, for flavor-neutral pions {i = j), there are disconnected 
contributions coming from one or more hairpin insertions. At LO, these appear only for taste 
singlet, vector, or axial-vector pions. Denoting the Minkowski hairpin vertices as we 
have jl^l : 

a^Sy, T= G {C,fj_} (taste vector); 

a^S'j^, Ts G {^/is} (taste axial- vector) ; ^ 

4mQ/3, T= = ^/ (taste singlet;) 

0, Ts G {^fiu,C,5} (taste tensor or pseudoscalar) 
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with 

16 
7 

The disconnected pion propagator is then 



K{A) = -J^{C2V{A) -C5V{A)) ■ (72) 



where 

^ i i jP^ - ^^,5) jP^ - (p^ - ^h) ^^^^ 

For concreteness we have assumed that there are three sea-quark flavors: m, rf, and s; the 
generalization to A^j flavors is immediate. Here = = is the mass of a taste-S 

pion made from a m and a u quark, neglecting hairpin mixing (and similarly for m£),H and 
^s,'e)i ^■n°,'Ei ^Ti,E, and m^/^s are the mass eigenvalues after mixing is included, and the ie 
terms have been left implicit. When specifying the particular member of a taste multiplet 
appearing in the disconnected propagator is unnecessary, we will abuse this notation slightly 
following Eq. and refer to 1^^^,^,, l^ui'i'^ '^in'i'- flavor-taste notation we have: 

I '^ab'^b'a' I {P) = \ ^ia,jl3^j'f^',i'a' I (p) = Sij6fi' V T~pTf ^^Vf^ ^^^, (75) 
I. J disc I. J disc ' ^ ' 

Equations (fTHjl . (|7Hl . and (fTHjl apply both to full and partially quenched QCD. In the 
latter case, i and i' can represent (quenched) valence quarks. Below, we will need the 
disconnected propagator only when i and i' both denote a particular valence quark: call it 
flavor X. Letting X be the flavor-neutral meson consisting of an x and an x, we will use the 
simplified notation Vxx = T^xx xx o'^) fo'^ example, I'xx = xx- 

We now calculate the one-loop decay constant, of a heavy-light pseudoscalar meson 
with light flavor-taste a. The calculation is done in the 4+4+4 theory, where there are four 
tastes per flavor. We then adjust the result for the 1 + 1 + 1 (one taste per flavor) case of 
interest. We work out explicitly the partially quenched case; full theory results are easily 
obtained by taking appropriate limits. The valence flavor is called x; the valence taste, a. 
In other words, we replace a ^ xa. As will be seen, the decay constant does not depend on 
the taste of the valence quark, so we may write /^^ Jb^- 

At one loop, we express the decay constant as 

/s. = (^1 + iQ^2p ^f^- + analytic terms^ . (76) 
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Recall that the lowest order term /^^ depends on the light valence flavor in a rather trivial 
way: /|^^ = /t / . The one- loop diagrams that contribute to SJb^ are shown in Figs. ^ 
and 121 where the cross indicates a disconnected pion propagator {i.e., one or more hairpin 
insertions for the singlet, vector and axial tastes). The terms in Fig. IHl vanish because the 
integrals produce an overall factor of on shell, which is then multiplied by the projector 
g^" - v^'v". 

The diagrams in Fig. contribute to the heavy-light wavefunction renormalization. Us- 
ing the propagators and vertices deflned above, one easily sees that the self energy S^b is 
proportional to the identity in taste space, and is diagonal in flavor. We thus have 

Sab{v-k) = Sia,jp{v-k) = SapSij Si{v-k) . (77) 

Flavor symmetry (broken by the diagonal mass terms) guarantees this flavor structure. The 
taste-independence follows from the residual discrete taste symmetry of the light quarks, 
which corresponds to the shift symmetry in the staggered action. This symmetry is 

gi (1 ® ^^,)qi , Qi qi{l ® ^^) , (78) 
at the level of the Symanzik action; while at the chiral level it is 



S 
a 
H 
H 



(79) 

Note that the symmetry is diagonal in flavor; the transformation acts only on the taste 
indices and affects all light quark flavors identically. 

The one-loop wavefunction renormalization, 6Zb^, for is then 

1 dSxiv-k) 



5Z 



v-k=0 



2 d{v-k) 

This contributes an amount Sf§^ to We flnd (in the 4+4+4 theory): 

1 



(80) 



9l 



E 



+ V'^^ + AVl^ + 4P 



XX 



(81) 
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where S runs over the 16 pion tastes, and / runs over the three sea quarks, u, d, and s. 
The term summed over / and S arises from the "connected" diagram. Fig. ^a); as we 
will see below, it involves an internal sea quark loop. The remaining terms come from the 
"disconnected" diagram. Fig. Q^b). We have summed over the four degenerate vector and 
axial disconnected contributions, as seen by the factors of four in front of and V^. 

The current corrections contribute an amount Sf^J to SJ'b^ ■ Computing the diagrams in 
Fig. 121 we obtain (in the 4+4+4 theory): 



cur 



1 



- + V' + \V\^ + 4P^ 



XX 



(82) 



These chiral integrals are similar to those that appear for the decay constant, /p, 
of a partially quenched pion P. The main difference is that here we have a single light 
valence quark x, whereas P is a bound state of two light valence quarks, x and y. So the 
corresponding integrals in the fp case have additional contributions with x ^ y la various 
terms. 

We now need to convert Eqs. ()81|) and ()82p from a 4+4+4 to a 1 + 1 + 1 theory. This can 
easily be done with a quark- flow analysis, following Ref. [isj. The g.,^ vertex at the quark 
level is shown in Fig. EJ^a); the current vertex at second order in the pion fields is depicted 
in Fig. ^b). The connected terms in Eqs. (j^?T|l and ()82j) that come from these vertices thus 
involve an internal quark loop, shown in Fig. Efa), and should be multiplied by a factor of 
1/4. As described in Ref. the only other changes in going from 4+4+4 to 1+1+1 appear 
in the mass eigenstates of the full fiavor-neutral propagators. This is due to the factors of 
1/4 that are associated with iteration of the two-point vertex, Eq. ()71|). in Fig. |^c). 

In the 1 + 1 + 1 theory, we thus have 



d^p 



p^p" 



[v ■ p + ie) 



4 Z-^ jp. 



It 



+ + \Vl^ + 4P^^ 



5/. 



cur 



d^p 



-T- 



+ + \V\^ + 41)^ 



XX 



(83) 



(84) 



where the mass eigenvalues that appear in T>^ ^ , and [c.f. Eq. ()74p ] are now the 1+1+1 
values For example, with three sea-quark fiavors, mrj/j ~ rriQ for large ttiq, rather than 
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the value 4mQ that follows from Eq. ()22|). After the ml oo limit is taken, this implies 
that the contribution is 4 times larger in the 1 + 1 + 1 theory than in the 4+4+4 theory. 

The disconnected propagators in Eqs. (jS^ and (jS^ can then be written as a sum of 
single or double poles using the residue functions introduced in Ref. We define {m} = 
{7711,7712, ... ,mn} as the set of masses that appear in the denominator of Eq. (|74j) . and 
{/i} = {/ii,/i2, . . . , fik} as the numerator set of masses. Then, for n > /c and all masses 
distinct, we have: 



J["'^l ({m};{/.}) 



nli(g^-/^f) 



E 



-1 



in+fc+l p[">fc] 



rrij + ze 



where 



771, 



({m};{/i}) = — 



If there is one double pole term for = m| (where G {^}), then 



(85) 



(86) 



^dp'^Hme; {m};{fi}) = 



(g2 - mj 

,2 



+ nj=i(^^ - + ie) 



9m« 



— 771 j + ze 



(87) 



We use the same definitions of the residues i?^"''^^ as in Ref. Ql; the factors of (— 1)"+*-'+! 
arise from the fact that the residues were defined for Euclidean space integrals, but we are 
working in Minkowski space here. 

Using Eqs. (jH^ and ()87|) . the integration in Eq. ()84p is immediate, and that in Eq. 



is standard in the continuum heavy-light literature. We follow Ref. 
dimensions, with chiral scale A^) 



3611 and define (in 4 — e 



J''^(m,A) = iAl 



{2tt)^-^ (p2 - m'^ i^^l^y . p - /\ 



-A[Ji(m, A)g^"' + J2(m, A^v' 



(88) 



A in this equation should not be confused with the A= defined in Eq. (jHEI)- For Eq. (jH^ . 
we need only 



dJf'^im, A) 



(9^ 



dA 

d dJ^"'{m,A) 



3 9AJi(m,A) 



dA 



A=0 



A=0 



167r2 



OA 



A=0 

3 d dJi{m,A) 



IGn"^ dTn? OA 



167r2 
3 



A=0 



167r 



£(m2) 
-Mm") 



(89) 
(90) 
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As in Ref. we use the arrow to denote the fact that we drop the analytic terms and 
keep only the chiral logarithms. The functions £ and ^ are defined as 

,2 



m In — — [infinite volume] 



d 



[m = 



dm? 



(m 



In 



m 



[infinite volume] 



(91) 
(92) 



These are the infinite volume forms; the effect of finite spatial volume will be discussed in 

Sec, ivrrn 

For Eq. ()84j] we need the Minkowski version of the integrals defined in Ref. [l^ : 

d^-'p 1 1 



him' 



him' 



iA 



^ J (27r)4-^~ (p2 _ m2 + ie) 



167r2 



[m 



167r2 



[m 



(93) 
(94) 



VI. FINAL NLO RESULTS 



Now we are ready to write down the full NLO result, including the analytic terms. 
After performing the integrals using Eqs. ()89|) through ()94|) . we get for the 1+1+1 partially 
quenched case with all masses unequal: 



fhO 
Ji 



1 l+H 1 



Bx / 1+1+1 



167r2/2 

: y 



-y 

16 ^ 



d 



dm\r L 



d 



dm\ 



i?fl(A^i.^'^);/ig))£(m^^) +[V^A] 



+ cJm„ + + mj + c^m^ + c^a 



(95) 



where / runs over the three sea quarks m, d, and s, S runs over the 16 meson tastes, and 
the sets of masses in the residues are defined as follows (with taste labels implicit): 



_y\^(3,x) 

_A/J(4,x) 



{m 



ml} 



r 2 2 2i 

{m^,m^o,m I 



{m 



;,mj,} 



(96) 
(97) 
(98) 
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The analytic terms in Eq. ()95p come from wave function renormalization contributions from 
£3 [Eqs. (Uni), and (j3Hj) ]. as well as from higher order corrections to the current ^'2'^ 
[Eqs. (HZj), (inD) and (jEH)]- The low energy constants (LECs) c^, and Ca can be written 
as linear combinations of the large number of constants appearing in £3 and However, 
such expressions are unlikely to be useful, and we omit them here. 

Of relevance to MILC simulations is the case where the up and down quark masses are 
degenerate, the 2 + 1 case. Setting m„ = rrici = mi, we obtain 



Bx 



fLO 

J Bx y 2+1 



1 + 



^91 



2 

d 



-E 

16 ^ 



3 ^ dm\ J 



— a 



d 



+ Cs{2mi + rris) + c^m^ + c^a^ 



(99) 



where 



{ml, ml} , 
{ml,ml} , 

r 2 2 2 1 

{mx,m^,m^,} 



(100) 
(101) 
(102) 



Another interesting limit is the case of Nf degenerate flavors of sea quarks of mass mj. 



Bx 



fLO 

■> Bx / N f degen 



1 + 



1 l + ^gl] Nf 



16 ^ 



[m 



xf,E) 



Nf 



{m\ j - mjf j) i{m\ j) - (.{m\j) 



d 



{m)f y - m\y) (i{m\ y) - {m)f y - ml, y) e{ml y) 



2 2 



+ [V ^ A]^^ + NfC^mj + c^m^ + c^a^ 



(103) 



where, in this case. 



m 



rj',V 



2 2r/ 



V ! 



(104) 
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and similarly for m^, ^. The result in Eq. ()lU3p reduces to that of Sharpe and Zhang 
Issj l in the continuum limit. ^ The 1/iV/ dependence of the taste-singlet term in Eq. p03|) is 
characteristic of the disconnected contribution and arises from the fact that rrijj'j ~ Nf m^/'i 
for large in the 1 + 1 + 1 + • ■ • theory. 

We now turn to full QCD, where the valence quark is one of the light sea quarks. In the 
2 + 1 case, we have 




(105) 



+ Cs(2m/ + nis) + c^rUs + CaO? . (106) 

We have used the fact that m^^ = {2m'^ j + m'^ j)/3 at this order to simplify the taste- 
singlet residues. Note that, as expected for full QCD, there are no longer any double-pole 
contributions. In the continuum limit Eqs. (jl()5|l and (jl()6j) reduce to known full QCD results 

0,0, Q. 

Note that the analytic term coming from the —1 in our Eq. ^ is absorbed into the LECs in Eq. (3.6) 
of Ref. 0- 
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VII. CHIRAL-SCALE DEPENDENCE AND LOW ENERGY CONSTANTS 



The effect of a change in the chiral scale can be absorbed into the LECs c^, c^, and 
Ca- In fact, the change of most of the one- loop terms is proportional to taste-split squared 
meson masses, rather than quark masses per se, since it is the meson masses that appear 
in the non-analytic terms. It has therefore proved convenient for light jl^ and heavy-light 
0] lattice calculations to rewrite the analytic terms as functions of meson masses, rather 
than quark masses, and, specifically, to use those combinations of meson masses that arise 
naturally from a change in scale in the non-analytic terms. Since squared meson masses have 
additive contributions proportional to due to taste violations [see Eq. (jUH|) ]. this procedure 
automatically introduces some terms proportional to Cy and Cg, and is equivalent to a 
redefinition Ca ^ Ca = Ca — PvCv — PsCs, where are constants. 

The redefinition is convenient because the scale dependence of Ca is simplified and no 
longer depends on the splittings As. More importantly, the procedure has a practical 
advantage for reducing lattice errors. From the simulations, Ca seems to be significantly 
smaller than Cq. In other words, most of the discretization error from the light quarks 
appears to be due to the dependence of the light meson masses in the chiral loops. To the 
extent this is true, i.e., that Ca can be neglected in first approximation, the dependence is 
tied to the mass dependence through and Cg. This means that a determination of the mass 
dependence at fixed lattice spacing can be used to estimate the lattice-spacing dependence. 

The appropriate redefinitions do not depend on the quark masses, but, like the LECs 
themselves, they do depend on Nf. Therefore it is most convenient to work with the case 
of degenerate sea quark masses but arbitrary Nj, Eq. (jl(J3|) . Using Eq. we can write 
the change in the chiral log contribution Sfs^ [c.f. Eq. (fTUj) ] under a change A^ in the 

chiral scale as 



6fB^ + log(AVA2) 



(l + 3^^) [ iV|-4 



) 





(107) 
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where 



with Aj the taste-singlet sphtting and A the average sphtting: 

A = 45"AH = 4(4^^ + 6Ar + 4Av + A,) . (110) 
lo lo 

Based on Eq. ()l()7j) . it is natural to redefine the LECs by 

NfCsTTif + Cymx + CaO? Nf Cs{2fimf + a^Asea) + C^(2/im^. + a^A^ai) + CaO^ (111) 

for degenerate quarks [e.g., in Eq. ()l()H|l ]. or 

c^(m„ + mrf + H ) + Cym^ + c^a^ {2fi {niu + + H ) + A^/a^Agca) 

+a(2/im^ + a^A^ai) + C^a^ (112) 

in the general case [e.g., in Eq. With these redefinitions, the change in the LECs 

under a chiral scale change is rather simple: 



,2 ^,2 ^1 + ^91) (^/-4) 



a(A,) = C.(AJ - log(A^/A^ 



327r2/2 



a(AJ = C.(AJ - log(A^/Aj) (5; + 5;,) . (113) 



We emphasize that the redistribution of the terms in Eq. pil|) or ()112j] is not unique, 
and it is possible to move them from Cg to without complicating the transformation 
properties of Ca under chiral scale change. In the "mixed" meson terms fn^f^-^ in Eq. fll03p . 
for example, we have chosen to associate half of each splitting with the valence quark x and 
half with the sea quark /. This is simple and "natural," but not required. The freedom 
stems from the fact that we have only one condition coming from chiral scale invariance, 
but two LECs {Cy and Cs) with which we want to associate a? dependence. This lack of 
uniqueness is particularly apparent in the case Nf = 2, where the choice we made for A^ai 
in Eq. ()1U8|) blows up. The problem is that, for Nf = 2, the change in the chiral logs 
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under change in chiral scale has no term that is proportional to rrix [see Eq. fjl07p ]. In 
this case, one might wish to associate all the splittings with the sea quark. Of course, a 
redefinition like Eq. or (|112|) is completely consistent for any choice of Agea and Aval- 

Some redefinitions are "better" than others, however, if they result in a smaller C^, and our 
choice does seem to work in the Nf = 3 case P|. 



VIII. FINITE VOLUME CORRECTIONS 



In computing the finite volume corrections, we assume that the lattice size in the time 
direction is larger than in the spatial directions, and indeed is large enough that it may 
3e treated as infinite to a very good approximation. This is the case for the MILC lattices 



IG 



25| . We are thus interested in corrections due to finite spatial volume only. Furthermore, 
we assume that the heavy-light meson is at rest: f'^ = (1,0,0,0). This is the case in all 
lattice simulations of the decay constant except those of "moving NRQCD" l22|. 

For integrals of the the form of Eqs. and (jM}, which come from light-meson tadpole 
loops, the corrections are expressed in a convenient form in Ref. jl^. The results of these 
integrals appear in our final answers, Eqs. through ()106|) . in the terms multiplied by 
the "1" in the factor 1 + 3g^. For such terms, we merely need to make the replacements 12 1 



(m^) £i(m^) = £(m^) + m^5i(mL) [finite volume, "1 terms"] , (114) 
(m^) £i(m^) = £(m^) + 53(mL) [finite volume, "1 terms"] , (US) 



in Eqs. (jHU and (HD), where 



^^(^L) = 7;t11 Ji > (116) 



5^{mL) = 2^ifo(|r|mL) , (117) 

with Kq and Ki the Bessel functions of imaginary argument. 

The integrals that involve a heavy-light meson in the loop appear in our final answers 
in the terms multiplied by g"^. For these integrals, Eqs. ()89j) and (jHOI), the finite volume 
corrections can be determined by comparison with the corrections to Eqs. ()93|) and ()94|) . To 
do this, we first perform the integrals over p° in both cases. This integration is not affected 
by finite volume corrections since the time extent of the lattices is taken to be infinite. We 
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obtain 



him') 



dA 



A=0 



(27r)3 2^p2 + m2 ' 
d^p p'^ 

(27r) 3 2(^2+^2)3/2 

Ji + 2m^ 



a^2_ • (118) 

Since the finite volume correction to Ji is given by Eq. ()114j) . the correction to the heavy- 
light integral Eq. ()89p is now determined via Eq. (jllSj) . Corrections to the integrals with 
double poles, Eqs. (jH^ and (jHOI), follow by differentiation with respect to rri'. 

The prescription for finite volume corrections to the terms proportional to g^, is then 

2 

i(m') ^ tgim') = i{m') + — (5i(mL) - 2b^{mL)) [finite volume, "^^ terms"] (119) 

3 

i{m') lg{n?) = i{m'') + 5^{mL) - ^55(mL) [finite volume, "^^ terms"], (120) 
where 

55(mL) = -—5'^{niL) = niL^ \r\Ki{\r\niL) , (121) 

with the prime denoting differentiation with respect to the argument. We have used Kq[z) = 
—Ki{z) and the relation 5^{mL) = —6i{mL) — {mL/2)5[{mL) [l^ . 

We note that our procedure for finding finite volume corrections of heavy-light integrals 
works well only for the rather simple integrals needed here. It does not easily generalize to 
cases where the time extent is finite or where the variable A in Eq. ()88|) is nonzero. The 
latter case is particular important since it is needed for semileptonic form factors, as well 
as for l/mq and 0{m'^) corrections. The finite size effects for heavy-light XPT have been 
studied in much more generality in Ref. [4^. 



IX. REMARKS AND CONCLUSIONS 



We have determined the heavy-light SXPT Lagrangian to NLO and used it to compute 
the heavy-light decay constant in partially quenched and full QCD. The heavy-light part of 
the chiral Lagrangian is identical at leading order to that in the continuum, aside from the 
extra taste degrees of freedom, which enter trivially. Thus taste-violating effects appear in 
the one-loop chiral logarithms only through the light-light meson sector, where the relevant 
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2^. 



parameters have already been determined by simulations [im2J,|25[. All the operators that 
appear at NLO in the heavy-light Lagrangian serve to produce only a term proportional to 

in the decay constants, and thus a single new low energy constant: Ca [Eqs. (jUHl) through 
flTOT)]) ] or Ca [Eq. (UnD or (fTT^ ]. 

We emphasize that our calculations are based on the assumption that the staggered 
"fourth root trick" is valid. We further assume that this trick is correctly implemented on 
the chiral theory by insertions of factors of 1/4 for each sea quark loop, where the loops are 
found by quark flow or replica analysis. 

Our results most relevant at present are those for the 2 + 1 case {Nf = 3 and mu = md), 
since they apply to calculations using existing MILC configurations for the light dynamical 
quarks. These results have proved quite useful in reducing the errors from the chiral and 
continuum extrapolations in a lattice calculation of the leptonic decay constants of D and 
Dg mesons j^. With the Lagrangian and the techniques presented here, one may also find 
expressions for the form factors for semileptonic heavy-to-light decays (such as D ^ Kiu 
or B ^ niu) in SXPT. This has been done jl^; a detailed report of the calculation is in 
preparation ML The result has already been used in Ref. 0]. We also call attention to the 
work of Laiho j4l|, who has applied heavy-light SXPT to study heavy-to- heavy semileptonic 
decays, e.g., B D*ih'. 

A calculation of heavy-light B parameters in SXPT should also be straightforward. Al- 
though additional operators will contribute to analytic terms at NLO, the one-loop chiral 
logs themselves will not involve any new low energy constants beyond those already present 
in the continuum or in the one-loop SXPT for light-light mesons. 

It is interesting to look at a simple example to see the effect of taste violations. In 
Fig. ini we plot fB^/fh'^ for the partially quenched case with two degenerate sea quarks: 
Eq. ()103p with Nf = 2. The solid line includes taste-breaking terms; while the dashed line 
shows the results with taste-breaking terms set to zero, i.e., the continuum limit. We choose 
parameters from the coarse MILC data set (a ~ 0.125fm). The light meson masses and 
splittings, as well as the values of 5^and 5^, come from measurements of the light hadron 
spectrum and decay constants jlO, 24, 25 1. We have used the value aruf = 0.010 for the sea 
quark mass. We set = = and choose to give a slope at large mass similar to what 
is seen in simulations. The continuum result shows the characteristic divergence of partially 



quenched chiral perturbation theory as the valence mass is taken to zero for fixed sea mass. 
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The taste violations, however, completely wash out this behavior and give a result that looks 
rather linear. This is in fact an example of an infrared sensitive quantity for which the chiral 



42] . Fitting to and removing the effects of 



and continuum limits in SXPT will not commute 
staggered taste violations is therefore crucial in controlling the systematic error associated 
with the extrapolation to the physical values of the light quark masses. 
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FIG. 1: Terms which contribute to the one- loop heavy-hght meson decay constant, arising from 
wavefunction renormahzation. The thick sohd hne is the heavy-hght meson, and the dashed hne is 
the hght pseudoscalar (i.e., "pion"). The sohd square is the current insertion. The cross indicates 
a disconnected propagator, Eq. ()75|): a pion hne without a cross is connected. 
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FIG. 2: Contributions to the one- loop heavy- light meson decay constant coming from corrections 
to the current insertion. Symbols have the same meaning as in Fig. ^ 
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FIG. 3: Terms which in principle arise at one loop in the decay constant but vanish since they are 
proportional to v'^(gi^i, — v^Vi,). 




FIG. 4: Relevant vertices at the quark level, (a) is the vertex, (b) is the second-order diagram 
coming from the current. 
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FIG. 5: Quark-flow diagrams that contribute to /b, corresponding to Figs. ^ (left hand side) and|2l 
(right hand side). The dashed lines here are the light quarks while solid lines are the heavy quarks. 
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FIG. 6: The ratio of the NLO decay constant for two degenerate sea quarks, given by Eq. H1U3|1 
with Nf = 2. The sohd hne is with the taste violations included, while the dashed line is the 
continuum limit. 
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